In-situ measurement of self-heating in intrinsic tunnelling spectroscopy. 
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Using advanced sample engineering we performed simultaneous measurement of interlayer tun- 
nelling characteristics and in-situ monitoring of temperature in Bi-2212 mesas. Together with a 
systematic study of size-dependence of interlayer tunnelling this allowed unambiguous discrimina- 
tion between artifacts of self-heating and gaps in electronic spectra of Bi-2212. Such a confident 
spectroscopic information, which is not affected by self-heating or surface deterioration, was ob- 
tained for the first time for a High-T c superconductor. We also derived general expressions and 
formulated main principles of self-heating valid for a large variety of materials. 
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PACS numbers: 74.50.+r 74.72.-h 74.25. Fy 44.10.+i 

Surface spectroscopy of high temperature super- 
conductors (HTSC) experience considerable difficulties 
caused by rapid chemical deterioration of the surface and 
a very short coherence length, due to which surface and 
bulk properties can differ at a scale of ~ one atomic layer. 
Those problems are avoided in intrinsic tunnelling spec- 
troscopy (ITS), which utilizes intrinsic Josephson junc- 
tions (IJJ's) naturally formed in highly anisotropic HTSC 
0. ITS has a superior resolution due to superconductor- 
insulator-superconductor structure of IJJ's and provides 
a unique opportunity to probe bulk electronic properties, 
indispensable for understanding HTSC |2j, La, |j, |5j . Un- 
fortunately, ITS is liable to self- heating caused by poor 
thermal conductivity of HTSC 0, 0, 0, E| • So far reports 
on self-heating in ITS differ by orders of magnitude even 
in similar short-pulse measurements 0, 0- In view of 
unique abilities of ITS, lack of consensus between differ- 
ent spectroscopic studies of HTSC, and recent contro- 
versy about the role of self-heating in ITS [ToL Hl| it is 
important to understand to what extent ITS is affected 
by self-heating. 

Self-heating in superconductors is practically unstud- 
ied, even though it is clear that it can be substantial 
due to inherently poor thermal conductivity at low T. 
For example, it was suspected for a long time that hys- 
teresis in Current- Voltage characteristics (IVC's) of over- 
damped Josephson junctions may indicate a substantial 
self-heating. However, it is still unknown to what ex- 
tent self-heating affects spectroscopy of superconductors 
even in more conventional cases of STM, photoemission 
or IVC's of low-T c junctions. Self-heating is also a grow- 
ing problem for semiconducting devices often having a 
mesa geometry [l2j. 

Here we present a comprehensive study of self-heating 
in Bi-2212 mesa structures. Advanced sample engineer- 
ing allowed simultaneous measurement of IVC's and T 
using a small portion of the mesa, nano-patterned by Fo- 
cused Ion Beam (FIB), as an in-situ thermometer. To- 
gether with a systematic study of size-dependence of ITS, 



this allowed unambiguous discrimination of self-heating 
artifacts from gaps in electronic spectra of Bi-2212. Such 
a confident spectroscopic information, not affected by 
self-heating or surface deterioration, was obtained for the 
first time for a HTSC. We also formulated main princi- 
ples and derived general expressions for self-heating in 
diffusive and ballistic cases. Both the developed exper- 
imental technique and analytic results can be used for 
studying self-heating in a large variety of materials. 

To distinguish artifacts of self-heating we first ana- 
lyze the shape of IVC's, assuming that electronic Density 
of States (DoS) is featureless, but the resistance exhibit 
strong dependence R(T) 



V = IR{T + AT). 



(1) 



Here T = T + AT is the mesa temperature, which is 
higher than the base temperature To due to self-heating. 
We consider a small circular mesa of radius a, containing 
N IJJ's, on top of a crystal of thickness t with thermal 
conductivities K a b,c in the ab— plane and c— axis, respec- 
tively. Heat flow has phononic and electronic channels. 
The latter is likely to be diffusive and allows the exact 
solution 0, which for a <C 2ty/ K a b/^ c is 
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where q is a power density. However, phononic transport 
is ballistic in the c-axis direction since phononic mean free 
path l p h ~ 1/xm |14| is much larger than the height of the 
mesa, containing few atomic layers. We estimate ballistic 
contribution assuming that phonons first shoot to the 
depth ~ l p h below the mesa and then spread via diffusion. 
In this case AT can be estimated from Eq.(3) of Ref. 6] 
with r = and z = l p h, multiplied by the fraction of 
heat, /, carried by phonons. The total self-heating is the 
sum of phononic and electronic contributions: 
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FIG. 1: Simulations for the self-heating model: a) dl /dV(V) 
and b) corresponding T(V) curves for four mesas with differ- 
ent geometry at To = 150^. A pronounced hump in dl/dV 
is seen for three largest mesas. Circles in Fig. b) indicate po- 
sitions of the humps. Irrespective of the sample geometry the 
hump occurs at the same T ~ 200K, eventhough the hump 
shape and voltage are different. Inset shows R(T) and k(T) 
dependencies used in simulations. 
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FIG. 2: Experimental characteristics at a) To = 16. 5K and 
b) To = 105K for mesas with different area fabricated on the 
same Bi-2212 single crystal, c) Size dependence of the peak 
voltage at To = 16. 5K and the hump voltage at To = 1057^. 
Dashed line shows the fit for ballistic heating . 



As expected, it coincides with Eq. for l p h — > 0. A 
large difference in anisotropics of electric p c /p a b ~ 10 5 
[3L fl3| and thermal K a b/ K c ~ 8 01 conductivities implies 
that c— axis heat transport in Bi-2212 is predominantly 
phononic, / ~ 1. Then we can write a simple estimation, 



AT bal 



qa 



4:K a bl 



ph 



(4) 



valid for a < l v h\J n a b/ n c . Here we added a factor 1/2 
since only ~ 1/2 of heat is going back to the mesa. 

To calculate IVC's we solved Eas. (|ll2[) with k = n(Tn + 
AT), using typical R(T) and K a b(T) dependencies [l4j, 
shown in inset of Fig. 1, and K a b/n c — 10. Fig. 1 shows 
simulated dI/dV(V) and T(V) curves for mesas with dif- 
ferent geometry. Humps in dI/dV(V) appear for three 
largest mesas. The origin of humps is easy to under- 
stand. At low To heating leads to a decrease of R. The 
ratio V/ 1 reaches minimum at the crossover temperature 
T* , see inset in Fig. 1. At higher bias, IVC's start to 
bend backwards because R(T) increases at T > T*. 

Even though the shape of self-heating IVC's depends 
on details of R(T) and k(T), sample geometry and heat 
transport mechanisms, artifacts of self-heating in all cases 
can be understood from the following three principles: 

1. T — q curves are single valued. At a given q mesas 
will be heated to higher T if started from higher Tq. 



2. For mesas with the same R{T), k(T) and T self- 
heating features in IVC's occur at the same T irrespective 
of mesa geometry. This is clearly seen in Fig. 1 b). 

3. Self-heating depends on the mesa size and the num- 
ber of junctions. From the second principle and Eas. l|2l4|l 
hump voltages in diffusive and ballistic cases are: 

Vdiff = V/N = Ki/Va~N cx N~ 1/2 A- 1/4 , (5) 
v b aii = V/N = K 2 /Va T N cx N^l 2 A^l 2 . (6) 

Here A is the mesa area and Kx^iTo) are material con- 
stants. The dependence Eq.© can be clearly traced from 
Fig. 1. This provides an unambiguous way to discrim- 
inate electronic spectra from self-heating artifacts: Self- 
heating depends on the sample geometry in contrast to 
the electronic spectrum, which is a material property. 

So far there was no systematic size-dependent ITS 
study. Since self-heating depends on R(T) and n(T), 
comparison should be made for mesas at the same crys- 
tal. Fig. 2 shows dI/dV(V) curves at T = 16.5K and 
105K for mesas with different area fabricated at the same 
optimally doped Bi-2212 single crystal, T c ~ 93K. All 
mesas contained = 9 IJJ's. Clearly distinguishable 
a peak at Tq < T c and a hump at Tq > T c are seen, 
previously attributed to the superconducting gap (SG) 
and the pseudo-gap (PC), respectively 0, §. Fig. 2 
c) shows peak, V p , and hump, Vh, voltages vs. mesa 
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FIG. 3: a) Thermometer mesa resistance as a function of 
the base temperature for several currents through the heater 
mesa. Inset shows the sample geometry, b) Normalized 
self-heating characteristics for several heating currents. In- 
set shows the IVC of the heater mesa at To = 5. 6 A". The 
actual temperature along the IVC is indicated. 



area. It is seen that for small mesas V p is independent 
of area, but for A > 16/xm 2 , V p becomes slightly smaller 
due to self-heating. This happens at P(V P ) > lmW. 
A stronger size-dependence is seen for the hump. For 
A > 10/im 2 , Vh is well described by the ballistic self- 
heating Eq.@), as indicated by the dashed line. However, 
for A < 10 fim 2 , Vh also shows a tendency for saturation. 
Here P(Vh) < 1.5mW. Note that scattering of data for 
small mesas is caused by flattening of the hump, see Fig. 
2 b). It is not clear if the residual flat hump represents 
the DoS, even though, STM data 0] and short pulse ex- 
periments |9(, which are less prone to self-heating, would 
imply so. Nevertheless, for such a flat hump it is more 
instructive to consider the characteristic slope of the V- 
shape dI/dV(V) at small bias, rather than Vh- 

From Fig. 2 it is seen that experimental ITS character- 
istics, retain the same V-shape and slope below the peak 
and the hump, irrespective of the mesa size. Such behav- 
ior is in stark contrast to self-heating artifacts shown in 
Fig. 1 a), for which both the slope and the voltage of the 
hump exhibited strong and correlated size-dependence. 
Size-independence of the shape of ITS characteristics to- 
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FIG. 4: Four-probe ITS characteristics at To slightly below 
and above T c ~ 93K. Sample geometry is shown in the top in- 
set. Bottom inset shows in-situ measured self-heating AT /P. 



gether with saturation of peak and hump voltages at 
A — > implies that those indeed represent gaps in DoS, 
rather than artifacts of self-heating. Such conclusion is 
consistent with reported strong suppression of the peak 
by magnetic field [5j, which is hard to explain in terms 
of self-heating, since k decreases with field 0. 

The ultimate judgement about the extent of self- 
heating can be made only by direct measurement of the 
mesa temperature. This requires fabrication of an in-situ 
thermometer situated in an intimate vicinity and having 
good thermal coupling to the mesa, and allowing inde- 
pendent calibration. So far there were no measurements 
which would satisfy all those requirements, even though 
several attempts has been made 0, El ■ 

To facilitate such measurements we developed a new 
method for fabrication of nano-scale multi-terminal ar- 
rays of IJJ's: first a mesa was made on top of a Bi- 
2212 crystal by self-alignment cross-bar photolithogra- 
phy. Then a small portion of the initial mesa was cut 
using FIB. The bias current was applied through the 
larger mesa (the heater) and the R(T) dependence of 
the nearby smaller mesa (the thermometer) was used for 
in-situ measurement of T Q. Typically we stopped the 
FIB cut within the height of the initial mesa so that both 
the heater and the thermometer mesas were sitting on a 
common pedestal - the bottom part of the initial mesa. 
The sketch of the sample geometry is shown in inset of 
Fig. 3 a) . Nano-scale separation between the heater and 
the thermometer mesas and small VT in the pedestal Q 
facilitated accurate measurement of self-heating. 

Fig. 3 a) shows R(Tq) of the thermometer for a set 
of dc-currents through the heater mesa, Iheat , for an un- 
derdoped Bi-2212 crystal, T c ~ 81K. It is seen that 
R decreases with Iheat as a result of heating. The AT 
can be obtained by subtracting To vs R curves at certain 
Iheat from that without heating, Iheat = 0. Correspond- 
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ing self-heating characteristics, normalized by the dissi- 
pation power are shown in Fig. 3 b). It is seen that at low 
Tq, AT/P decreases with Iheat- This is due to a strong 
k(T) dependence at low T, see inset in Fig. 1. At higher 
To > 40 — 50 A, AT/P curves for different Iheat approach 
each other, indicating flattening of k(T) dependence. 

Inset in Fig. 3 b) shows the IVC of the heater mesa at 
T = 5.6A'. The actual T along the IVC is indicated. It 
is seen that at the end of sum-gap knee (above the peak) 
the mesa is heated by 19. 2A, i.e. considerably smaller 
than T c . The peak voltage is not affected by such heating 
since the SG has a flat dependence in this T— range [J 

Sample geometry used in this work is also interest- 
ing because it facilitates true four-probe measurements 
of IJJ's in the pedestal, removing questions about the 
influence of contact resistance between the mesa and the 
electrode and improving ITS resolution. Fig. 4 shows 
four-probe characteristics from 90A to 302A. Here we 
can clearly observe the peak just few A below T c ~ 93A 
as well as some, presumably strong-coupling phonon fea- 
tures above the SG, which disappear at T c . Inset in Fig. 4 
shows in-situ measured self-heating AT/P for this mesa, 
which is similar to that in Fig. 3 b) and all other mea- 
sured samples. In Fig. 4 it is indicated that self-heating 
at the peak at To = 90 A is small AT ~ 0.5 A even despite 
a large number of IJJ's, N = 90 in total. Thus, influence 
of self-heating on the peak vanishes at To — » T c because 
both AT/P decreases with T and the power at the peak 
vanishes together with the SG at T — > T c |j| . The realis- 
tic simulation of how self-heating affects T— dependence 
of the SG can be found in Ref.y|. 

From Fig. 4 it is seen that dI/dV(V) are V-shaped at 
T > T c with a pronounced dip at V = 0. Whether or not 
the dip is a self-heating artifact can be understood from 
in-situ T— measurements, shown in inset. In Fig. 4 it is 
indicated that at Tq = 90A the mesa is heated to ~ I00A 
at V = 200ml/. On the other hand, at T = 100 A the 
self- heating free dI/dV(V — 0) ~ 4.0mS is considerably 
smaller than dI/dV(V = 200mV) ~ 14.8mS. Therefore, 
self-heating is insufficient to explain such a large increase 
of conductance from V = to 200mV and the V-shape 
characteristics must be attributed to the PG in DoS. 

Now we can reanalyze one of the most important re- 
sults of previous ITS studies - evidence for coexistence 
of the SG and the PG 0, 0, 0, which favors different 
origins of the two gaps. When the peak and the hump 
are observed simultaneously, the hump occurs at larger 
V and, therefore, at higher T than the peak. From Fig. 
2 it is seen that for large mesas the hump is affected by 
sclf-hcating and is exaggerated due to progressive back- 
bending of dI/dV(V) at large V. In most of previous ITS 
experiments dealing with larger mesas, humps were cer- 
tainly affected by self-heating. This is particularly true 
for Bi-220I compound ^H, which has a similar AT/P 
as Bi-2212, but considerably lower T c and, thus, more 
prone to self-heating. So, does the PG exists at T < T c ? 



The answer can be obtained from Fig. 4. It is seen that 
the dl/dViy) curve at T = 90 A acquires characteris- 
tic V-shape right after the peak, where self-heating is still 
small. Thus, the PG does coexist with superconductivity 
at T < T c . Such conclusion is also supported by observa- 
tion of V-shape PG characteristics at T < T, in magnetic 
fields H > H C 2 5] and in the vortex core fig . 

In summary, we performed a comprehensive analysis 
of self-heating in Bi-2212 mesas, which allowed unam- 
biguous discrimination of gaps in electronic spectra from 
artifacts of self-heating. This was achieved via systematic 
size-dependent study and in-situ measurement of self- 
heating using a nano-patterned part of the mesa as the 
in-situ thermometer. We observed that for small mesas 
both shapes and voltages of the peak and the hump (or 
rather V-shaped suppression of conductance at V = 0) 
are independent of the mesa size and appear at T consid- 
erably smaller than T c and T*, respectively. Therefore, 
they are not due to self-heating 10], but represent the 
superconducting gap and the normal state pseudo-gap 
in electronic DoS, respectively. Such a confident spec- 
troscopic information for a HTSC not affected by self- 
heating, or surface deterioration was obtained here for 
the first time. On the other hand, ITS characteristics 
of larger mesas can be strongly affected by sclf-hcating. 
The reported threshold mesa size and dissipation power, 
at which self-hcating becomes insignificant, as well as the 
measured values of self-heating AT/P are typical for our 
mesas, but not universal. They depend on the sample 
geometry, materials and experimental setup. Thus, it is 
important to carefully design samples for ITS: decrease 
mesa sizes, number of junctions, avoid suspended struc- 
tures and, in particular, to employ the top heat spreading 
layer, which is the most important heat sinking channel 
in mesas 0, ^2 • 

Finally we want to note that derived expressions and 
formulated principles of self-heating are general and valid 
for any material. Similarly, the developed experimental 
method, in which a small portion of the sample is used 
as the in-situ detector, can be used for analysis of self- 
heating in a large variety of materials. The only require- 
ment is that the resistance of material or the contact 
between electrode and material should have considerable 
T— dependence, while layered structure and HTSC is not 
essential. For example, our method can be directly ap- 
plied for studying self-heating in small mesa-like semi- 
conducting transistors [l2j |. 
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ERRATUM 

The units of the vertical axis in Fig. 3a) should read 
10fcf2, not fcO. This does not affect any part of the 
manuscript because only relative shapes of the curves 
were used for determining shifts in temperature. 



